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Abstract. Starting from the Larmor clock we introduce a hierarchy of density of 
states. At the bottom of this hierarchy are the partial densities of states which rep- 
resent the contribution to the local density of states if both the incident and the 
out-going scattering channel is prescribed. We discuss the role of the partial densities 
of states in a number of electrical conduction problems in phase coherent mesoscopic 
systems: The partial densities of states play a prominent role in measurements with 
a scanning tunneling microscope on multiprobe conductors in the presence of current 
flow. The partial densities of states determine the degree of dephasing generated by a 
weakly coupled voltage probe. We show that the partial densities of states determine 
the frequency-dependent response of mesoscopic conductors in the presence of slowly 
oscillating voltages applied to the contacts of the sample. We introduce the off-diagonal 
elements of the partial density of states matrix to describe fluctuation processes. These 
examples demonstrate that the analysis of the Larmor clock has a wide range of appli- 
cations. 



1 Introduction 

The Larmor clock is one of the most widely discussed approaches to determine 
the time-scales of tunneling processes. The essential idea[^]j^,^ of the Larmor 
clock is that the motion of the spin polarization in a narrow region of magnetic 
field can be exploited to provide information on the time carriers spend in this 
region. It is assumed that incident carriers are spin polarized and that they im- 
pinge on a region to which a small magnetic field is applied perpendicular to the 
spin polarization of the incident carriers (see Fig. |). The spin polarization of 
the transmitted and reflected carriers can then be compared with the polariza- 
tion of the incident carriers. Dividing the angle between the polarization of the 
exciting carriers and that of the incident carrier by the Larmor frequency ujl 
gives a time. Originally, only spin precession (the movement of the polarization 
in the plane perpendicular to the magnetic field) was considered. However, Ref. 

pointed out, that especially if we deal with regions in which only evanescent 
waves exist (tunneling problems) the polarization executes not only a precession 
but also a rotation into the direction of the magnetic field. In fact in the presence 
of a tunneling barrier, the spin rotation, is the dominant effect. Ref. Q consid- 
ered a rectangular barrier and considered a magnetic field of the same spatial 
extend as the barrier. In the local version of the Larmor clock, introduced by 
Leavens and AersQ, we consider an arbitrary region in which the magnetic field 
is non- vanishing and investigate again the direction of the spin polarization and 
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rotation of the transmitted and reflected carriers. The magnetic field might be 
non-vanishing in a small region localized inside the barrier or in a small region 
outside barrier on the side on which carriers are incident or on the far side of 
the barrier. We mention already here, that the response of the carriers is highly 
non-local: Even carriers which are reflected are affected by a magnetic field that 
is non- vanishing only on the far side of the barrier where naively we would expect 
only transmitted carriers B. In this work we use the local Larmor clock to derive 




Fig. 1. Spin polarized carriers incident on a barrier subject to a weak magnetic field 
Bo- The transmitted carriers exhibit both a spin precession and and a spin rotation. 
After g . 

a set of local density of states states |^J^,§,§[@ which we call partial densities 
of states which are related to spin precession and in terms of sensitivities which 
are related to spin rotation. The partial densities of states, below abbreviated as 
PDOS, are useful to understand a number of transport problems: the transmis- 
sion probability from a tunneling microscope tip into a multiterminal mesoscopic 
conductor can be expressed in terms of PDOS, the absorption of carriers by 
an optical potential (a potential with a small imaginary component), inelastic 
scattering and dephasing caused by a weak coupling voltage probe, and the low 
frequency transport in mesoscopic conductors. 

The partial densities of states are determined by functional derivatives of 
the scattering matrix P,[7p|,p^ . Only in certain limited situations can the den- 



Larmor clock, Density of States 



3 



sity of states be expressed in terms of energy derivatives. Expressions for the 
density of states in terms of energy derivatives of the scattering matrix are fa- 
mihar In the discussion of characteristic times the distinction between 

time-scales found from energy derivatives (like the Wigner-Smith phase delay) 
and time scales found from derivatives with respect to the local potential (the 
dwell time) has found some recognition. In contrast density of states are almost 
invariable discussed in terms of energy derivatives. Here we emphasize that a 
more precise discussion of density of states also uses derivatives with respect to 
the (local) potential and not energy derivatives. It is the dwell time (or sums of 
relevant dwell times) which are related to the density of states |13;1J|. The use 
of energy derivatives always signals that approximations are involved. 

The interpretation of the Larmor clock remains a subject of discussion. Ref. 
[D considered the total rotation angle dived by the Larmor precession frequency 
to be the relevant time. This interpretation brings the Larmor clock into agree- 
ment with the time-scales obtained by considering tunneling through a barrier 



with an oscillating potential |15 . Subsequent works have argued that the pre- 
cession angle and the rotation angle dived by the Larmor precession frequency 
separately should be viewed as time scales The difficulty with such an 

interpretation is not only that one has two scales characterizing the same pro- 
cess, but the times defined in such a way are also not necessarily positive. Since 
we aim at characterizing a duration, that is a definite draw back. The two time 
scales can be combined into a complex time, with the real part referring to the 
precession time and the imaginary part to the time obtained from rotation. Like 
negative times, complex times are not part of the commonly accepted notions of 
time. Steinberg argues that the clock presents only a "weak measurement" and 



that therefore complex times are permitted 1 18 . In quantum mechanics the ques- 



tions "how much time has the transmitted particle spent in a given interval" is 
problematic since being in an "interval" and "to be transmitted" corresponds to 
noncommuting operators p9[ | . Reasonably, we can only speak of a time duration, 
if it is real and positive. 

A comparison of the Larmor clock with the closely related linear ac-response 
of an electrical conductor shows immediately the ambiguity of the clock: in the 
ac-response of a conductor which is predominantly capacitive (tunneling limit) 
the voltage leads the current whereas for a highly transmissive conductor the 
response is inductive and the current leads the voltage. 

The language used here implies similarly an extension of the usual notion of 
density of states. At the bottom of the hierarchy of density of states which we 
discuss are the partial densities of states (PDOS) which represents the contri- 
bution to the local density of states if we prescribe both the incident and the 
out-going channel. It turns out that certain partial densities of states are not pos- 
itive. (These are of course just the PDOS that correspond to negative precession 
times). Thus the discussion presented below does not solve the interpretational 
questions related to the Larmor clock. Nevertheless as we will show even such 
negative PDOS are physically relevant. Using the partial densities of states, ei- 
ther by summing over the out-going channels (or by summing over the incident 
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channels) we obtain the injectivity of a contact into a point within the sample or 
the emissivity of a point within the sample into a contact. Both injectivities and 
emissivities are positive and in the language of tunneling times correspond to 
local dwell times for which either the incident channel (or the outgoing channel 
but not both) are prescribed. Finally, if we take the sum of all the injectivities 
or the sum of all the emissivities we obtain the local density of states. 

2 The scattering matrix and the local Larmor clock 

We start by considering a one-dimensional scattering problem [^|[|^. We con- 
sider particles moving along the y-axis in a potential V{y). The potential is 
arbitrary, except that asymptotically, for large negative and large positive val- 
ues of y it is assumed to be flat. We adopt here the language from mesoscopic 
transport discussions and call the region of large negative y the contact 1 (the 
left contact) and the region of large positive y contact 2 (the right contact). We 
assume that the quantum mechanical evolution is described by the Schrodinger 
equation. A particle with energy E has for large negative y a wave vector ki {E) 
and a velocity vi (E) and for large positive values of y a wave vector ^2 (E) and 
a velocity V2{E). We are interested in scattering states. A particle incident from 
the left is for large negative values of y described by a scattering state 

^i(£;,y) = e^'=^^-f ^ne-^^-^^ (1) 

and for large positive values of y is described by a transmitted wave 

ME.y)^{-f'^S2ie^^y. (2) 

Similarly, a particle incident from the right is for large positive values of y de- 
scribed by a scattering state 

^2(i?,y) = e-^'"^ + ^22e-^'=^^ (3) 
and for large negative values of y is described by a transmitted wave 

ME.v)^{-Y'^S,2e^'''y. (4) 

Vl 

Here the amplitudes Safi determine the elements of the 2 * 2-scattcring matrix 
of the problem. Each scattering matrix element is a function of the energy E 
of the incident carrier and is a functional of the potential V{y). To express 
this dependence we write SapiE, V{y)). Conservation of the probability current 
requires that this matrix is unitary and in the absence of a magnetic field time- 
reversal invariance implies that it is also symmetric. 

Next we now consider a weak magnetic field applied to a small region. The 
magnetic field shall point into the z-direction. For simplicity we consider the 
case where the magnetic field is constant in a small interval [y, y + dy] and 
takes there the value B. We consider only the effect of the Zeeman energy. 
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Thus the motion of the particle remains one-dimensional and is as in the ab- 
sence of a magnetic field confined to the y-axis. In the set-up for the Lar- 
mor clock, we consider particles with a spin. The spin of the incident particles 
(in the asymptotic regions) is polarized along the x-axis. For spin 1/2 parti- 
cles the wave functions are now spinors with two components ip^{y,E) and 
tp-{y,E). Carriers incident from the left, have a spinor with components given 
by ip+{y,E) = 'tp-{y,E) = {l/^/2)exp{ik\y). The Zeeman energy which is gen- 
erated by the local magnetic field is diagonal in the spin up and spin down 
components. Consequently, in the interval [y, y + dy] for a particle with spin up, 
the energy is reduced by Tiojl/'2' with uji,j2 = giiB/% and for a particle with spin 
down the potential is increased by %ljl/2. Thus with the magnetic field switched 
on particles with spin up travel in a potential V{y) — dV{y) and particles with 
spin down travel in potential V{y) + dV{y). Here V{y) is the potential in the 
absence of the magnetic field and dV{y) is the potential generated by the mag- 
netic field. Thus dV{y) vanishes everywhere, except in the interval [y,y + dy] 
where it takes the value dV{y) = %u)l/'^ = gfJ-B. 

We can evaluate the polarization of the transmitted particles and the reflected 
particles if we can determine the scattering matrix for spin up and spin down 
particles in the potential generated by the magnetic field. Thus we need the 
scattering matrices S^p{E, V{y) =F dV{y)) where 5*+ is the scattering matrix for 
spin up carriers and S" is the scattering matrix for spin down carriers. Since 
the potential variation generated by the magnetic field is small we can expand 
these matrices away from the scattering matrix for the unperturbed potential. 
Thus we find to first order in the magnetic field for the scattering matrices 



S^^{E, V{y)) T dV{y)) = S^p{E, V{y)) T [dSo^0{E, V {y)) / dV {y)\dV {y)dy + ...(5) 



The variation of the scattering matrix due to the magnetic field is proportional 
to the derivative of the scattering matrix with respect to the local potential at 
the location where the magnetic field is non-vanishing. More generally, we can 
consider a magnetic field which varies along the y-axis (but always points along 
the z-axis). This leads to a potential 5V{y) determined by the local magnetic 
field. The variation of the scattering matrix is then determined by a functional 
derivative [5sa[j{E,V{y))/5V{y)\ of the scattering matrix with regard to the 
local potential. 



We emphasize that even though this equation looks quite simple, the evaluation 
of a functional derivative of a scattering matrix, while not difficult for a one- 
dimensional problem, can still be a laborious calculation. 

Let us now find the precession and rotation angles of the polarization of 
the transmitted and reflected carriers. The normalized spinor of the transmitted 
particles which determines the spin orientation of the transmitted a particles has 



S%{E,V{y))T6V{y)) = S^p 



{E,V{y)) 
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the components 

First consider the polarization in the y direction. It is found by evaluating the 
expectation value of the Pauli spin matrix (jy, 



< S >2^^ ^ < ^J^ >= _j^^l^2^^1^. (9) 

^ 2 2 (1^2+12 + 1^2- |2) ^ ' 



Here the indices 21 indicate that we consider transmission from left (1) to right 
(2) and evaluate the the spin in the transmitted beam. We need the spin po- 
larization only to first order in the applied magnetic field. Using Eq. we 
find 

< sy >2i= -iy{2, y, l)ujLdy (10) 

where T — |>S'2iP is the transmission probability in the absence of the magnetic 
field and 

is the partial density of states at y of carriers which emanate from contact 1 (the 
asymptotic region for large negative y) and eventually in the future reach contact 
2. Since initially the spin polarization was in the x-direction < Sy >2i directly 
determines the angle of precession of the carriers in the x — y-plane. Thus by 
dividing < Sy >2i by the Larmor precession frequency we can formally introduce 
a quantity with the dimension of time, which we call Tj,(2, y, 1) and which is given 
by Ty{2, y, 1) — [h/T)v{2, y, l)dy. Here the index y indicates that we deal with a 
time-scale obtained from the y-component of the spin polarization. We can now 
proceed to evaluate also the y-component of the spin polarization of the carriers 
which are refiected and can proceed to evaluate the y-component of the spin 
polarization of the carriers that arc in the past incident from contact 2 (large 
positive y) and in the future will be transmitted into contact 1 (large negative 
1) or will be reflected back into contact 1. We can summarize the results in the 
following manner: There are a total of four spin polarizations to be considered, 
each of them determined by a partial density of states 



.(a,y,/?) = --^ Si,^ - -^S^0 (12) 




4ni y °''^6V{y) 

of carriers that are incident at contact /3 = 1, 2 and eventually in the future are 
transmitted or reflected into contact a — 1,2. Formally, the time scales related 
to precession in the local magnetic field at y can be introduced which are related 
to the partial densities of states via, Ty{a,y,l3) — {h/\Sai3\'^)i'{a,y,(3)dy where 
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I'S'q./sP is the transmission probability T if a and /3 are not equal and is the 
reflection probability i? if a and f3 are equal. Thus with each element of the 
scattering matrix we can associate a partial density of states. Later we discuss 
the properties of the partial densities of states in more detail. 

Next we consider the spin polarization in the z-direction. For the carriers 
incident in contact 1 and transmitted into contact 2 we find that the z-component 
of the transmitted carriers is determined by 

>2i=-<^ikJVi>=-^^t^— (13) 



Using Eq. (||) we find 



where we call 



<s,>i=^i]{2,y,l)LULdy. (14) 



the sensitivity of the scattering problem. Since the spin polarization of the in- 
cident particles was originally along the x-direction only a small z-coniponent 
of the incident particle determines a spin rotation angle. We can formally in- 
troduce a time scale Tz{2,y,l) associated with spin rotation which is given by 
Tz{2, y, 1) = {h/T)ri{2, y, l)dy. Again we can ask about the z-polarization of re- 
flected particles and can ask about the z-polarization of particles incident from 
the right. The results are summarized by attributing each scattering matrix 
element Safj a sensitivity 

•"°--« = -i(^i»i| + ^*-) 

which determine the time-scales Tz{a,y,(3) — {h/\Sa(3\'^)r]{a,y, P)dy. Finally we 
can determine the spin polarization in the x-direction. This component is reduced 
from its initial value both because of spin precession in the x ~ y-plane and 
because of the rotation of spins into the z-direction. Since we have at every 
space point 

<Sx>^ + < Sy >2 + < .s^ >2= n:^/4:, (17) 

it follows immediately that the time scale Tx is related to the two time-scales 
introduced above by 

rx^ir'y+T^z)'/'. (18) 

Using the expressions for Ty and Tz given above, we find for the time-scales 
the following expressions 
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We reemphasizc that neither the partial densities of states nor the sensitivities 
are in general positive. In contrast, Tx{a^y, (3) is positive for all elements of the 
scattering matrix. 

3 Absorption and Emission of Particles: Injectivities and 
Emissivities 

Before discussing the partial densities of states in more detail it is of interest to 
investigate the absorption of particles in a small scattering region j2^. We as- 
sume that in a narrow interval [y, y + dy] there exists a non- vanishing absorption 
rate F . Thus the potential V{y) is equal to Vb(2/) — ifiF in the interval [j/, y + dy] 
and is equal to Vo{y) out-side this interval. To solve the scattering problem we 
need to find the scattering matrix S^^{E, V{y)) in the presence of this complex 
potential V{y). Of interest here is, as in Ref. the limit of small absorption. 
The case of strong absorption is also of interest but thus has been used only to 
discuss global properties and not the local quantities of interest here |^,|3) . For 
a small absorption rate we can expand the scattering matrix S^p{E,V{y)) in 
powers of the absorption rate away from the scattering problem in the original 
real potential Vq. We obtain 

+in[5Sap{E, V{y))/SV{y)]\viy)^v„(y)rdy + (20) 

We note that the adjoint scattering matrix has to be evaluated in the potential 
V*{y) and hence 

S^J<iE,V{y)) ^ Si^{E,Vo{y)) 

-zn[6Sl^{E, V{y))/SV{y)]\y^y^=v,^y-^rdy + .... (21) 

With these results it easy to show that the transmission and reflection probabil- 
ities in the presence of a small absorption in the interval [y^y + dy] are 

\S!:p{E,V{y))\' = \S^f,iE,Voiy))\^il-r,.{a,y,(3)dy) (22) 

where {SapiE, Vb(j/))p is the transmission probability T of the scattering prob- 
lem without absorption if a and /3 are different and is the reflection probability 
R of the scattering problem without absorption if a = /3. The incident current 
jin, must be equal to the sum of the transmitted current Jt, the reflected current 
jn and the absorbed current jr-, 

jin = jr+jB+jr- (23) 

Using Eq. (^0|) and taking into account that the incident flux is normalized to 
1, we find for carriers incident from the left or right /3 = 1,2 an absorbed flux 
given by 

jriy,f3)^riyiy,(3)dy, (24) 
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where v[y^f3) is called the injectivity of contact /3 into point y. The injectivity 
of the contact is related to the partial densities of states via 

u{y,P) = J2Ha.y,P)- (25) 

a 

In our problem with two contacts the injectivity is just the sum of two partial 
densities of states. 

Another way of determining the absorbed flux proceeds as follows. The ab- 
sorbed flux is proportional to the integrated density of particles in the region 
of absorption (in the interval [y,y + dy]). The density of particles can be found 
from the scattering state "0/3 (y) given by Eqs. ^ -^). For carriers incident from 
contact (3 the absorbed flux is thus 

jr{y,(3)=rj^^\My)\'dy. (26) 

Note that here the density of states 1/hvp of the asymptotic scattering region 
appears. It normalizes the incident current to 1. Thus we have found a wave 
function representation for the injectivity. Comparing Eq. (^) and Eq. ( p2|) 
gives 

K2/,/3)= (27) 

The total local density of states vly) at point y is obtained by considering carriers 
incident from both contacts. In terms of wave functions v{y) is for our one- 
dimensional problem given by 

Ky)-ET^I^'3(y)l'- (28) 

Thus the total density of states is also the sum of the injectivities from the left 
and right contacts 

13 

There is now an interesting additional problem to be addressed. Instead of a 
potential which acts as a carrier sink (as an absorber) we can ask about a 
potential which acts as a carrier source. Obviously, all we have to do to turn our 
potential into a carrier source is to change the sign of the imaginary part of the 
potential. With a a carrier source in the interval [y, y + dy] we should observe 
a particle current toward contact 1 and a particle current toward contact 2. We 
suppose that carriers are incident both from the left and the right and evaluate 
the currents in the contact regions. The total current injected into the sample 
at y is 

jin{y) ^ rv{y)dy. (30) 

Taking into account that the incident current is normalized to 1 the current 
joutiP, y) in contact [3 due to a carrier source at y is given by 

jout{a,y) = \~Y.K3{E,V{y))\'' (31) 
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due to the modification of both the transmission and reflection coefhcients. Using 
Eq. ( pO[ ) (with r replaced by -F) gives 

jout{a,y) = -r^iy{a,y,f3)dy = -rv{a,y)dy. (32) 

The current in contact a is determined by the emissivity i^i^a, y) of the point y 
into contact a. Note the reversal of the sequence of arguments in the emissivity 
as compared to the injectivity. Thus the emissivity is like the injectivity a sum 
of two partial densities of states, 

'^ioi,y) ^^i'{a,y,P). (33) 

For a scattering problem in the absence of a magnetic field the injectivity and 
emissivity are identical. If there is a homogeneous magnetic field present they 
are related by reciprocity: the injectivity from contact a into point y is equal to 
the emissivity of point y into the contact a in a magnetic field that has been 
reversed, j/+_B(y,a) = v^B{.a,y). 

We have thus obtained a hierarchy of density of states: At the bottom are the 
partial densities of states v{a, y, 0) for which we describe both the contact from 
which the carriers are incident and the contact through which the carriers have 
to exit. On the next higher level are the injectivities v{y^a) and the emissivities 
v{a,y). For the injectivity we prescribe the contact through which the carrier 
enters but the final contact is not prescribed. In the emissivity we prescribe 
the contact through which the carrier leaves but the incident contact is not 
prescribed. Finally, on the highest level is the local density of states i^(y) for 
which we prescribe neither the incident contact nor the contact through which 
carriers leaves. 

For simple scattering problems (delta- functions, barriers) the interested reader 
can find a derivation and discussion of partial densities of states in Refs. [p| jic|p0|| . 

Returning to time scales: we have shown that the partial densities of states are 
associated with spin precession. It is tempting, therefore, to associate them with 
a time duration. However, as can be shown, the partial densities of states are not 
necessarily positive. (The simple example of a resonant double barrier shows that 
one of the diagonal elements ^{a, y, a) has a range of energies where it is negative 
[pT|). The injectivities and emissivities are, however, always positive. The proof is 
given by Eq. (p5|). We can associate a local dwell time r£i(y, /?) = h,h'n{y, (i) with 
the injectivity which gives the time a carrier incident from contact (3 spends in 
the interval [y, y-\-dy] irrespective of whether it is finally reflected or transmitted. 
Similarly, we can associate a dwell time with the emissivity which is the time 
carriers spend in the interval [y, y + dy] irrespective from which contact they 
entered the scattering region. There is little question that the dwell times have 
the properties which we associate with the duration of a process: they are real 
and positive. However, as explained they do not characterize transmission or 
reflection processes. 
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4 Potential Perturbations 

Thus far our discussion has focused much on the partial densities of states. The 
sensitivity introduced as a measure of the spin rotation in the Larmor clock has, 
however, also an immediate direct interpretation. We have seen that the partial 
densities of states are obtained in response to a complex perturbation of the orig- 
inal potential V{y). The sensitivity comes into play if we consider a reaZ perturba- 
tion SV of the original potential. Thus if we consider a potential which is equal to 
V{y) + SV in the interval [y, y + dy] and equal to V{y) elsewhere the transmission 
probability in the presence of the perturbation is = T+Anri{a, y, (3)5Vdy, 
with Of ^ /3. The reflection probability is = R + 4TTri{a,y,a)SVdy. Since also 
+ BX = 1 we must have ri{y) = ri{a,y,/3) = —ri{a,y,a). For our scattering 
problem, described by a 2 * 2 scattering matrix there exists only one indepen- 
dent sensitivity r]{y). In the Larmor clock the sensitivity corresponds to spin 
rotation and the fact that there is only one sensitivity follows from the con- 
servation of angular momentum: the weak magnetic field cannot produce a net 
angular momentum. If carriers in the transmitted beam acquire a polarization in 
the direction of the magnetic field then carriers in the reflected beam must have 
a corresponding polarization opposite to the direction of the magnetic field. In 
mesoscopic physics, in electrical transport problems, the sensitivity plays a role 
in the discussion of non-linear current- voltage characteristics and plays a role if 
we ask about the change of the conductance in response to the variation of a 
gate voltage |Q. Below, we will not further discuss the sensitivity, but we will 
present a number of examples in which the partial densities of states play a role. 



5 Generalized Bardeen formulae 

It is well known that with a scanning tunneling microscope (STM) we can mea- 
sure the local density of states . STM measurements are typically performed 
in a two terminal geometry, in which the tip of the microscope represents one 
contact and the sample provides another contact . Here we consider a differ- 
ent geometry. We are interested in the transmission probability from an STM 
tip into the contact of a sample with two or more contacts as shown in Fig. 
Thus we deal with a multiterminal transmission problem If we denote the 
contacts of the sample by a Greek letter a = 1,2,.. and use tip to label the con- 
tact of the STM tip, we are interested in the tunneling probabilities TaUp from 
the tip into contact a of the sample. In this case the STM tip acts as carrier 
source. Similarly we ask about the transmission probability Tupa from a sample 
contact to the tip. In this case the STM tip acts as a carrier sink. Earlier work 
has addressed this problem either with the help of scattering matrices, electron 
wave dividers, or by applying the Fermi Golden Rule. Recently, Gramespacher 
and the author JlOf have returned to this problem and have derived expressions 
for these transmission probabilities from the scattering matrix of the full prob- 
lem (sample plus tip) . For a tunneling contact with a density of states i^up which 
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couples locally at the point x with a coupling energy \t\ these authors found 

Tup.a ^ 4,7^^ mp\t\'^ v{x,a) , (34) 

Tatip = 4TT^v{a,x)\t\'^iytip- (35) 

In a multiterminal sample the transmission probability from a contact a to 
the STM tip is given by the injectivity of contact a into the point x and the 
transmission probability from the tip to the contact a is given by the emissivity 
of the point x into contact a. Eqs. (^) and (|3^) when multiplied by the unit 
of conductance e'^/h are generalized Bardecn conductances for tunneling into 
multiprobe conductors. Since the local density of states of the tip is an even 
function of magnetic field and since the injectivity and emissivity are related by 
reciprocity we also have the reciprocity relation Ttip,a{B) — Ta^upi— B) . 

The presence of the tip also affects transmission and reflection at the massive 
contacts of the sample. To first order in the coupling energy |ip these probabil- 
ities are given by 

\S'7p\^ = M-47T'Hc^,x,P)\t\'n.p. (36) 

The correction to the transmission probabilities a ^ (3 and reflection probabili- 
ties a = /3 is determined by the partial densities of states, the coupling energy 
and the density of states in the tip. Note that if these probabilities are placed in 
a matrix then each row and each column of this matrix adds up to the number 
of quantum channels in the contacts. 




Fig. 2. Partial density of states measurement with a scanning tunneling microscope. 
The tip of an STM couples at a point x with a coupling strength t to the surface 
of a multi-terminal conductor. The contacts of the conductor are held at potentials 
/ia = eVa and the tip at potential /is — eVup. After [|l(|. 



6 Voltage probe and inelastic scattering 

Consider a two probe conductor much smaller than any inelastic or phase break- 
ing length. Carrier transport through such a structure can then be said to be 
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coherent and its conductance is at zero temperature given by G = {e^ /h)T, 
where T is the probabiUty for transmission form one contact to the other. How 
is this result affected by events which break the phase or by events which cause 
inelastic scattering? To investigate this question Ref. proposes to use an ad- 
ditional (third) contact to the sample. The third probe acts as a voltage probe 
which has its potential adjusted in such a way that there is no net current flowing 
into this additional probe, /s — 0. The current at the third probe is set to zero 
by floating the voltage /X3 = eV^ at this contact to a value for which vanishes. 
The third probe acts, therefore, like a voltage probe. Even though the total cur- 
rent at the voltage probe vanishes individual carriers can enter this probe if they 



are at the same time replaced by carriers emanating from the probe |26 . En- 
tering and leaving a contact are irreversible processes, since there is no definite 
phase relationship between a carrier that enters the contact and a carrier that 
leaves a contact. In a three probe conductor, the relationship between currents 
and voltages is given by la = '^pGapVp where the Gaf} are the conductance 
coefficients. Using the condition /a = to find the potential V-i and eliminating 
this potential in the equation for I2 or Ii gives for the two probe conductance 
in the presence of the voltage probe 

G = -G.,--^Hi^. (37) 

For a very weakly coupled voltage probe we can use Eqs. (^- ^6|). Taking into 
account that Gap — —(e^//i)| 5*0/3 p for a 7^ /3 we find 

G = ^ (t - 4.^ |tP [K2, 1) - ^^^^^^1^]) . (38) 

Here v{x) is the local density of states at the location of the point at which 
the voltage probe couples to the conductor. Eq. ( |3^ ) has a simple interpretation 
[ p6| . The first term T is the transmission probability of the conductor in the 
absence of the voltage probe. The first term inside the brackets proportional to 
the local partial density of states gives the reduction of coherent transmission 
due to the presence of the voltage probe. The second term in the brackets is 
the incoherent contribution to transport due to inelastic scattering induced by 
the voltage probe. It is proportional to the injectivity of contact 1 at point x. A 
fraction v{2^x)/v{x) of the carriers which reach this point, proportional to its 
emissivity, are scattered forward and, therefore, contribute to transport. Notice 
the different signs of these two contributions. The effect of inelastic scattering 
(or dephasing) can either enhance transport or diminish transport, depending 
on whether the reduction of coherent transmission (first term) or the increase 
due to incoherent transmission (second term) dominates. 

Instead of a voltage probe, we can also use an optical potential to simulate 
inelastic scattering or dephasing. However, in order to preserve current, we must 
use both an absorbing optical potential (to take carriers out) and an emitting 
optical potential (to reinsert carriers). The absorbed and re-emitted current must 
again exactly balance each other. From Eq. (p3) it is seen that the coherent 
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current is again diminished by ri>(2, x, 1), i. e. by the partial density of states at 
point X. The total absorbed current is proportional to Fv^x, 1), the injectance 
of contact 1 into this point. As shown in section 3 a carrier emitting optical 
potential at x generates a current — x) in contact 1 and generates a current 
—rv{2,x) in contact 2. It produces thus a total current —rv{x). In order that 
the generated and the absorbed current are equal we have to normalize the 
emitting optical potential such that it generates a total current proportional to 
rv{x^ 1) (equal to the absorbed current). The current at contact 2 generated by 
an optical potential normalized in such a way is thus —rv(2, x)v{x^ l)/v{x). The 
sum of the two contributions, the absorbed current and the re-emitted current 
gives an overall transmission (or conductance) which is given by Eq. ( ^8|) with 
47r|tp replaced by F. 

Thus the weakly coupled voltage probe (which has current conservation built 
in) and a discussion based on optical potentials coupled with a current conserving 
re-insertion of carriers are equivalent |^ . There are discussions in the literature 
which invoke optical potentials but do not re-insert carriers. Obviously, such 
discussions violate current conservation. A recent discussion pst , which compares 
the voltage probe model and the approach via optical potentials, does re-insert 
carriers but does this in an ad hoc manner. In fact Ref. claims that the 
Onsager symmetry relations are violated in the optical potential approach. This 
is an incorrect conclusion arising from the arbitrary manner in which carriers 
are re-inserted. 

We conclude this section with a cautionary remark: We have found here 
that the weakly coupled probe voltage probe model and the optical potential 
model are equivalent. But this equivalence rests on a particular description of 
the voltage probe. There are many different models and even in the weak coupling 
limit our description of the voltage probe given here is not unique. The claim can 
only be that for sufficiently weak optical absorption and re-insertion of carriers 
there exits one voltage probe model which gives the same answer. Differing weak 
coupling voltage probes are discussed in Refs. p9|l . 



7 AC Conductance of mesoscopic conductors 

In this section we discuss as an additional application of partial densities of 
states briefly the ac-conductance of mesoscopic systems. We consider a conduc- 
tor with an arbitrary number of contacts labeled by a Greek index a = 1, 2, 3.... 
The problem is to find the relationship between the currents /q(w) at frequency 
uj measured at the contacts of the sample in response to a sinusoidal voltage 
with amplitude Va(aj) applied to contact fi. The relationship between currents 
and voltages is given by a dynamical conductance matrix Gapito) such that 
-fa(w) = '^i3Gafi{i^)Vp{Lo). All clcctric fields are localized in space. The overall 
charge on the conductor is conserved. Consequently, current is also conserved 
and the currents depend only on voltage differences. Current conservation im- 
plies GajS = for each (3. In order that only voltage differences matter, the 
dynamical conductance matrix has to obey Gaf3 — for each a. We are 
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interested here in the low frequency behavior of the conductance and therefore 
we can expand the conductance in powers of the frequency j^], 

G^p{uo) = G% - iujE^p + K^pLu^ + 0(^3). (39) 

Here G^^ is the dc-conductance matrix. Eap is called the emittance matrix and 
governs the displacement currents. gives the response to second order in 
the frequency. All matrices G°^, Ea/3 and Kaf3 are real. 

We focus here on the emittance matrix Eafj. The conservation of the total 
charge can only be achieved by considering the long-range Coulomb interac- 
tion. Here we describe the long-range Coulomb interaction in a random phase 
approach (RPA) in terms of an effective interaction. The effective interaction 
potential g{x' , x) has to be found by solving a Poisson equation with a non-local 
screening term. The effective interaction gives the potential variation at point x' 
in response to a variation of the charge at point x. With the help of the effective 
interaction we find for the omittance matrix ||6|l 

Eap = 

Here the first term, proportional to the integrated partial density of states, 
is the ac-response at low frequencies which we would have in the absence of 
interactions. The second term has the following simple interpretation: an ac- 
voltage applied to contact f3 would (in the absence of interactions) lead to a 
charge built up at point x given by the injectivity of contact /3. Due to interaction, 
this charge generates at point x' a variation in the local potential which then 
induces a current in contact a proportional to the emissivity of this point into 
contact a. The effective interaction has the property that at an additional charge 
with a distribution proportional to the local density of states gives rise to a 
spatially uniform potential, / dx'v(x')g{x' ,x) — 1 for every x. This property 
ensures that the elements of each row and each column of the emittance matrix 
add up to zero. In particular if screening is local (over a length scale of a Thomas 
Fermi wave length) we have g{x',x) = 5{x' — x)v'^^{x). In this limit the close 
connection between Eq. ( ^ ) and Eq. (^5|) is then obvious. 



dxv{a,x,f3)— I dx' dxv{a,x') g{x' ,x) ^{x, (3) 



(40) 



8 Transition from Capacitive to Inductive Response 

The following example |^o|] provides an instructive application of the ac-conduc- 
tance formula Eq. (^). Consider the transmission through a narrow opening 
shown in Fig. |^. Carrier motion is in two dimension through a potential U(x, y) 
which has the form a saddle with a height Uq. The conductance (transmission) 
through such a narrow opening (a quantum point contact) has been widely 
studied and is found to rise step- like |3l| as a function of the potential Uq with 
plateaus at values G = {2e'^/h)N corresponding to perfect transmission of N spin 
degenerate channels. Here we are interested in the emittance £^ as a function of 
Uq. We consider the case that Uq is so large that transmission is completely 
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Fig. 3. Charge dipole across a saddle point constriction, /xi = eVi and fi2 = eV2 are 
the potentials ol the contacts, 5Ui and 5U2 axe local potentials. The dashed lines are 
the equipotential lines of the equilibrium potential. After [po|. 



blocked and then lower Uq such that the probabihty of transmission probabiHty 
T gradually increases from to 1. 

We introduce two regions Sli and [^2 to the left and the right of the barrier, 
respectively. Instead of the local partial density of states we consider the par- 
tial density of states integrated over the respective volumes i7i and f22- Thus 
we introduce Dakp = Jq^ dxdyv(a,x,y, (3). We furthermore introduce the total 
density of states D of the two regions. We assume that the potential has left- 
right symmetry and consequently the density of states in the regions fii and i?2 
are Di = D2 = D/2. Ref. |^ evaluates the partial densities of states semiclassi- 
cally. We find that carriers incident from contact 1 and transmitted into contact 
2 give rise to a partial densities of states in region 1 given by I?2ii = TDi/2. To 
determine D212, we note that in the semiclassical limit considered here, there 
are no states in f2i associated with scattering from contact 2 back to contact 2, 
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hence it holds D212 = 0. With similar arguments one finds for the semi-classical 
PDOS 

Dc^kp = Dk (T/2 + 6^p{R Sck - T/2) ) , if a, /3 ^ 3 . (41) 

From Eq. ( |4l| ) we obtain for the emissivity into contact 1 from region 1 and injec- 
tivity from contact 1 into region 1, Df^ — D\i — (1/4)(1 + R)D and and obtain 
for the emissivity into contact 1 from region 2 and the injcctivity into region 
1 from contact 2, = ^^12 — {^/^)TD. Instead of the full Poisson equation 
the effective interaction |3^,^,^ is determined with the help of a geometrical 
capacitance C. For a detailed discussion we have to refer the reader to Refs. 
[ ^0|j3^ . Due to charge conservation we have E = En = E22 = —E12 = —E21 
with 

E = {RC^ - (42) 

Here C^^ — R^^{C^^ + {e'^D/A)^^) is the effective capacitance of the contact. It 
is proportional to the reflection probability and proportional to the series capac- 
itance of the geometrical capacitance C and the "quantum capacitance" e^D. 
If the contact is completely closed we have R = 1,T = and the emittance is 
completely determined by the capacitance C^. If the channel is completely open 
we have R = 0,T = 1 and the emittance is E = —D/4. It is negative indicating 
that for a completely open channel the ac-response is now not capacitive but 
inductive. Thus there is a voltage Uq for which the emittance vanishes. For a 
simple saddle point potential the behavior of the capacitance and emittance is 
illustrated in Fig. ^. The dotted line shows the conductance, the dashed line is 
the capacitance and the full line is the emittance as a function of the saddle 
point potential Uq- The emittance is capacitive (positive) for a nearly closed 
contact and changes sign as the transmission probability increases from near 
to 1. The emittance shows additional structure associated with the successive 
opening of further quantum channels. 

A similar transition from capacitive to inductive behavior is found in the 
emittance of a mesoscopic wire: Guo and coworkers investigate the emit- 
tance of a wire as a function of impurity concentration. A ballistic wire with no 
impurities has an inductive response, a disordered metallic diffusive wire has a 
capacitive response. Experiments on ac-transport in mesoscopic structures are 
challenging and we can mention here only the work by Pieper and Price |]3^ on 
the ac-conductance of an Aharonov-Bohm ring and the recent work by Desrat 
et al. on the low-frequency impedance of quantized Hall conductors. 

Our simple example demonstrates the difficulty in associating a time with a 
result obtained by analyzing a stationary (or as here) a quasi-stationary scatter- 
ing problem. The emittance divided by the conductance quantum e'^/h has the 
dimension of a time. (In the non-interacting limit e^/C << 4/Z? we have pGf 
E/ (e^//i) = (i? — T)td/4: where tjj is the dwell time in the two regions f2i and 
/22). But since this "time" changes sign such an interpretation is not appropri- 
ate. Furthermore, in electrical problems the natural "times" are i?C-times if the 
low frequency dynamics is capacitive or an i?/L-time if it is inductive. The fact 
that we describe here a crossover from capacitive to an inductive- like behavior 
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demonstrates that neither of these two time-scales can adequately describe the 
dynamics. 




Fig. 4. Conductance of the saddle point constriction (dotted line) in units of 2e^//i, 
capacitance and emittance (dashed and full curve) as a function of the height of the 
saddle point potential C/q. In the range of voltages shown three quantum channels are 
opened. After Isol. 
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9 Partial Density of States Matrix 

Thus far the main aim of our discussion has been to iUustrate the concept of 
partial densities of states with a number of simple examples. We now would like 
to point to some important extensions of this concept. 

Quantum mechanics is a theory of probability amplitudes. We can thus sup- 
pose that our initial state is a scattering experiment which is described by a 
superposition !?■(?/) — ap^fj{y) + a^^^{y). Here and are two scattering 
states describing particles incident in channel /3 and in channel 7 with ampli- 
tudes afj and a-y. Often such superpositions are eliminated since we can as- 
sume that each incident amplitude carries its own phase (ji/j and (p^. We sup- 
pose that a/3 = \af3\^/'^exp{i(t)i3) and = |a^|"'^/^ea;p(i(/)^) and that these phases 
are random and uncorrelated. Thus if an average over many scattering exper- 
iments is taken we have < e*"^''e~*'^T >= Sp^y and consequently we find that 
< |!f'(y)P >= |a/3p|iZ'/3p -I- |a^p|lf'^p is the sum of scattering events in the dif- 
ferent quantum channels. However, as soon as we are interested in quantities 
which depend to forth (or higher order) on the amplitudes then even after av- 
eraging over the random phases of the incident waves wc find that scattering 
processes which involve two (or more) incident waves matter. For instance con- 
sider < |!f'(j/)P|!f'(2/')P > - < W{y)? >< W{y')? > which is determined by 
|a/3p|a^ptp'^(y)!f'^(y)tp'^(y')if'*(?/') -I- /i.e.. This expression describes a correlation 
of the particle density at the points y and y' . Obviously for such higher order cor- 
relations superpositions are very important even if we associate a random phase 
with each incident scattering states and take an average over these phases. 

We can relate such density fluctuations to functional derivatives of scattering 
matrix expressions. To describe density fluctuations of scattering processes with 
an incident carrier stream from channel /3 and 7 we consider 

A,). -3L(s:„^-^S„,) (43) 

Note that two different scattering matrices enter into this expression. The con- 
nection to the scattering states is determined by the relation 

^i.(a,j/,/3,7) = il/h)ivpv^)-'/'W*piy)i'^iy) (44) 

a 

Here Vfj and Vj are the (asymptotic) velocities of the carriers in the scattering 
channels /? and 7 far away from the scattering region. Eq. ( |44[) w as given in Ref. 
\ ^7\ and a detailed derivation of this relation is presented in ||38| . 

The expressions i/(a,y,/3,7) can be viewed as the off-diagonal elements of 
a partial density of states matrix. In Eq. (^) we take the sum over out-going 
channels. The resulting matrix is the local density of states matrix. Using Eq. 
(HI) in Eq. wc find, 

Hy, 7) - E A 7) = - ^ E {'IwS)) (45) 
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where we have taken into account that the scattering matrix is unitary. 

Let us now consider the total density of states matrix. The fluctuations of 
interest are then the total particle number fluctuations in the scattering region 
Q, D{(3,^) — ^ ^-^dyvijj, (3,^). Furthermore, if the volume of integration is suffi- 
ciently large, we can, in WKB approximation, replace the functional derivative 
with respect to V with a derivative with respect to energy. The matrix which 
governs the fluctuations in the particle number in the scattering region then 
becomes, 

Eq. ( ^ ) is the Wigner-Smith delay time matrix We have earlier emphasized 
that the appearance of the energy derivative is a consequence of approximations 
(here the fact that we consider the WKB limit). We also mention that strictly 
speaking, here we do not consider a "delay" . We do not compare with a reference 
scattering problem (a free motion) as is typically done in nuclear scattering 
problems. Eq. (^) determines a total time or absolute time and we should more 
appropriately call it the absolute time matrix instead of the delay time matrix. 

The partial density of states matrix has been used in Ref. [p7| to obtain the 
second order in frequency term of the ac-conductance (see Eq. (|39|)). Ref. po| 
investigated the current induced into a nearby gate due to charge fluctuations 
in quantum point contacts and chaotic cavities. More recently, the charge fluc- 
tuations in two nearby mesoscopic conductors was treated with this approach 
and the effect of quantum dephasing due to charge fluctuations was calculated 
within this approach ||4l|,|38| . The results can be compared with other theoretical 



works and with experiments |43 



The Wigner-Smith delay time matrix has received wider attention. In recent 
years the focus has been on the calculation of the entire distribution function 
of delay times jiij for structures whose dynamic is in the classical limit chaotic 
(chaotic cavities). Predominantly structures have been investigated in which 
carrier propagation is an allowed energy region. 

We conclude by briefly discussing the Wigner-Smith matrix for a tunnel 
barrier. For a symmetric barrier with transmission and reflection probability 
T and R the scattering matrix has elements Sn = S22 = —iVRexp{i4i)^ and 
S21 = S12 = y/T exp{i4>) where 4> is the phase accumulated during a reflection or 
transmission process. Thus the elements of the Wigner-Smith delay time matrix 
Eq. (H) are 

„ „ 1 d(bn „ „ 11 dTn , , 

I'll = 2?22 = T^i2 = T^-^i = : — -. (47) 

To be speciflc consider now the case of a tunnel barrier. In the WKB limit we 
have i? « 1 and T = exp{-2S/h) with 5" = j dy^/2m^JV{y) - E where the 
integral extends from one turning point to the other. We have d(j)/dE = and 
using the expression tt = m J dy{^/2m\JV{y) — E)~^ for the traversal time of 
tunneling gives dT/dE = {2TT/h)T. Consequently, the diagonal elements of the 
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Wigner-Smith matrix vanish and the non-diagonal elements are 

2?i2 = V21 = ^VT. (48) 
Znn 

Thus while the the average density inside the barrier vanishes (the trace of the 
Wigner-Smith matrix is zero) the off-diagonal elements are non-zero and indicate 
that the fluctuations of the charge will in general be non-vanishing even deep 
inside the classically forbidden region. 

The discussion of this section rests admittedly vague. While important first 
steps have been made to extend the notion of partial density of states to treat 
fluctuations, our discussion shows that even at the conceptual level, there is 
clearly room for more research. 

Another development which could be discussed here is a theory of quantum 
pumping in small systems. In quantum pumping one is interested in the current 
generated as two parameters (like gate voltages, magnetic fluxes) which modulate 
the system are varyed sinusoidally but out of phase. Brouwer p5[ , Avron et al. 

Shutenko et al. and Polianski and Brouwer |Q develop a theory which 
is based on the modulation of the partial densities of states discussed here. 

10 Discussion 

The Larmor clock and its close relatives have become one of the most widely 
investigated approaches mainly in order to understand the question: "How long 
does a particle traveling through a classically forbidden region (a tunnel barrier) 
interact with this region?" . We have already in the previous sections pointed 
out that there is no consensus in the interpretation even of this simple clock. 
Regardless of these difficulties the investigation of the Larmor clock has been 
helpful in understanding a number of transport problems: In particular we have 
discussed a hierarchy of density of states as they occur in open multiprobe meso- 
scopic conductors. These density of states are directly related to local Larmor 
times. We have shown that a small absorption or a small emission of particles 
can be described with these densities (or in terms of the Larmor times) . We have 
shown that the transmission probabilities through weakly coupled contacts like 
the STM is related to these densities. We have shown that a weakly coupled volt- 
age probe, describing inelastic scattering or a dephasing process can be treated 
in terms of these densities. We have also pointed out that the ac-conductance 
of a mesoscopic conductor at small frequencies can be expressed with the help 
of these densities. Furthermore, we have indicated that it is useful to consider 
also the off-diagonal elements of a partial density of states matrix since this 
permits a description of fluctuation processes. Thus there is no question that 
the investigation of the Larmor clock has been a very fruitful and important 
enterprise. 
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